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On the Part of the Parallactic Inequalities in the Moon's 

Motion which is a Function of the Mean 

Motions of the Sun and Moon. 

By Eenest W. Brown, M. A., Haverford College, Pa. 



In Vol. I of the American Journal of Mathematics,* Mr. Gr. W. Hill, taking 
into account only inequalities depending on the mean motions of the sun and 
moon, has shown that these inequalities are able to be determined tp a high 
degree of accuracy by using moving rectangular axes, one of which passes 
through the mean place of the sun. This paper is an adaptation of his method 
so as to include that class of inequalities which depends also on the ratio of the 
solar and lunar distances, and in particular, the principal part of the Parallactic 
Inequality. Owing to the use which has been made of this latter in obtaining 
the parallax of the sun, it becomes of importance that its coefficient should be 
accurately known. Delaunay's expressionsf are deficient in this respect, and we 
have no means of knowing how near Hansen's numerical value is to the truth. 
The inequalities are obtained below in an algebraical and numerical form, the 
latter giving their coefficients in longitude and parallax true to about one thou- 
sandth of a second of arc. 

I. 

Transformation of tlie Equations of Motion. 

The inequalities from purely circular motion which depend only on the 
lunar eccentricity and mean motions of the sun and moon only, are given by the 
equations 



dtf + Zn dt + r* y ~°' 



(1) 



* Researches in the Lunar Theory, pp. 5, 129, 245. t Mem. Fr. Acad. So., Vol. XXIX, p. 847. 
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where n', (i have their usual significations, the axes being rectangular and revolv- 
ing with uniform angular velocity n', and that of x passing through the mean 
sun. The earth is supposed to be moving in a circle round the sun with uniform 
angular velocity n'. The above equation is that which would be obtained from 
the consideration of a disturbing body of infinitely great mass, and at an infinite 
distance moving round the earth in a circle with finite uniform angular velocity 
n', where n'* = mass -i- (dist.) 3 . In order to get the inequalities free from the 
lunar eccentricity, Mr. Hill gets a particular solution of the above equations. 
When now we wish to include inequalities dependent on the distance of the sun, 
we can no longer suppose it at an infinite distance and of infinite mass. We must 
include the part of the disturbing function due to this cause which has been 
omitted. It will then be seen that we can obtain all the inequalities which a 
disturbing body moving round the earth with constant angular velocity in the 
moon's orbit produces, when the moon's undisturbed orbit is taken to be circular. 
Hence we have, instead of zero on the right-hand sides of equations (1), to 

put -j-1 and -t-^ respectively, where 

7% P 1 5 1 

+ ~aJ r [f ~ 5x$ y* + ~8~ xy J -*" 

fiiQT (Mil 

Multiplying the equations thus formed by —rr , -jr , adding and integrating 
the result, we get the Jacobian equation 

(£)'+ (£)-£-*"= «*-«'• 

Transforming to new coordinates u, s, where 

« = «; + yV — 1, 
s = x — yV — 1 , 
our equations of motion and the Jacobian integral become respectively 

_ F + afl V-i.- ar + f £ 5I i,=- Tfl r(u + .) + a- 3 ^: 

^i.^-^^ + ^iu + sf +2^-20, 
dt dt (us)* 4 v ' ' 1 ' 



(2) 
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and fix takes the form 

ni= IT ["16" ^ + sS ) + ~W m ( u + s )] 

+ 

n'W 3 
I2 X , expressed in terms of u, s, is obtained directly by expanding , where p 

P 

is the distance between the sun and moon and is equal in this case to 



V(a' a — 2a!r cob 6 + i*). 

For since r 2 = «s and 2r cos 6 = 2x = m + s , we have 2r n cos »0 = u n + s n , and 

therefore nV , r r * s 

-— = n'a'{l + lf P l + ir P,+ ....). 

P n is the zonal harmonic of degree n, and is expressible in terms of 

cos n6 , cos (n — 2) 6 ... . 



Hence we have 



where 



«'' . „ w' a 



1 a' 6 a' * 



_ 1.3....(2n-l) 1.3....(2n-3) J_ w _ 2 , . , 

rjP "~ 2.4.... 2n (M +*;+ 2.4. ...(2n— 2) 2 W ^ +S ) + •••• 

II. 

Solution of the Equations. 

In order to solve our equations for the particular class of inequalities which 
we wish to get, take the particular integrals 

x = 2 J.j cos iv (t — t ) , 
y = 2A { smiv(t — t ), 

i having positive integral values from zero to infinity, and v = n — »'. Put 

■A-i = <%i— i + a_ i_ i, B i -^a t _x — «_{_i, 
we obtain x = 2a 4 _ x cos iv (t — t ) , 

y = 2a i _ l sin iv (t — t ) , 
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where the summation is now extended to negative values of i. Transforming to 
the complex variables u, s, and putting e" ( ' - *» 11 ' _1 = £, we have 

u = Xa t _i^, 
Also putting 



m -v> x -rf'^dZ- t7~" '~dt~' 

the equations (2) may be written 

[*• + 2mD - !=?] «=- 4 ->+•>- » ! •-£■• f 5 ' 

[»-»*-$]> = --§-rf<.+.)—P. Jr-^, 

(tts)" 4 v ' n' 

Multiply the first of these equations by s, the second by u, and add to the last; 
also with the same multipliers subtract the first from the second, the resulting 
equations will be 

Z> 2 (us) — DuDs — 2m (uDs — sDu) + -J- m 3 ( M + s) s 



4 

2 / (ZQj rffli 



D {uDs - sDu - 2mus) + -|- m 3 (« 2 -s 2 ) = - m a . ~ («^ -s^ 1 ) ' 



(3) 



The constant * has vanished and will have to be determined from our previous 
equations in terms of the new constants introduced in the Particular Solution 
(Hill, p. 132). 

Substituting in (3) the values of u, s, which are 

u = X&i _ jf, s = l li a_ { _ g , 

or, what is the same thing, 

so that 

us = XiXj-ajaj _ £*, uDs — sZ>m = 2*2,- (i — 2/ — 2) afy _ £ ' , etc., 
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and equating the coefficient of £* to zero (except when i— 0), we obtain 

3, \f - U + 1)(» —J— 1) - 2m (t - 2; - 2) + -|- m 3 ] ofl,_, 

9 



+ -j- m% [a^xai^.i + flL-,_ 1 a_ <+> _J = — m 3 A-, 



(4) 



2,- [t (» — 2/ — 2) — 2mi] a/^_ 



3 



where — m 2 i ( , — m 3 !^ are the coefficients of £* on the right-hand sides of equa- 
tions (3). Multiply these equations by 2 and 3 respectively, and take the sum 
and difference, we get 

Xj [5i 2 - 8 (J + 1) i + 2 (/ + I) 2 - 2m (5» - 4j- 4) + 9m 3 ] aj a_ i+j 

+ 9m%o / _ 1 a i _,_ 1 = — m 2 (2£, + 3if f ) , 

2, [— » 3 + 4 (i + 1) i + 2 (/ + l) 3 + 2m (» -f 4y + 4) + 9m 3 ] a^_ <+ ,, 

+ 9m 3 2^_ j _ 1 a_i + ,._ 1 = — m 3 (2L, — 3M t ) . 

The terms of lowest order in these equations are a a t and a a_ t . In order 
to separate out these, multiply the last pair of equations by 

— ^ 3 + 4*+2+2m(* + 4) + 9m 2 , 
_ 5J2 + 8» — 2 + 2m (5* — 4) — 9m 3 , 

respectively, add the products and divide the whole by 

12* 2 [2 (i 2 — 1) — 4m + m 3 ] . 

In the resulting equation the term a a_ t will have vanished and the coefficient 
of a a t will be — 1 . The equation is 

2j { [*>/] a fl-i+i + H aj-iOt-j-i + (*) «-,_ift-«-w-it j 

where r . .-, _ JL (^--2)y + ^+ 2^-2-2 (y-i +2) m + m 3 

L*» ^ J i • 2(i 3 — 1) — 4m + m 3 

3m 3 i 3 — 4t — 2 — 2(t + 4) m — 9m 3 
W— — "4?" 2(i 2 — 1) — 4m + m 2 ' 

__ 3m 2 5J 2 — 8t + 2 — 2 (5» — 4) m + 9m 2 
0) — — 4^ 2(t 3 — 1) — 4m + m 3 ' 

15 
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and therefore, as they should be, 

[t, %] = — 1 [i, 0] = 0. 

The equation (5) corresponds to that obtained by Mr. Hill on p. 135 of his 
memoir referred to. It may be shown that if we put zero instead of the func- 
tion on the right-hand side of equation (5), every a with an odd suffix will 
vanish, and with one or two changes in notation it becomes the same as his. 

The expressions L t , M t must now be obtained. Since Q. x is formed of homo- 
geneous functions of u , s of the 3 a , 4 th degrees, we have 

d^ d£i r n" r 5 „ 3 -i 

+ 6 • "^ [256 ( w5 + s5) + 0& us ^ + s ^ + jk wV (« + s )] + • • • • ' 

dClx dil x n'° r 15 . 8 .. , 3 . -1 



u 



and -g- Zj, — K~M t being the coefficients in these when for u, s are substituted 
their values. Hence 



9 



- -i- [^ (u 3 ), + 4 (s 3 ), + B t (tAr)i + 5 4 ' K)J 

- i [4 (A + <?*' (**)i + A W + D'i (™%+ E i (tiW),] 



i- [**« (« B )< + ^ (s 6 ) 4 + 0, (tt**), + G[ {m% + H t (tA») ( + HI (t*V) J , 
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where (m 3 )*, (s 3 ) 4 , etc., denote the coefficients of £' in u 8 , s s , etc. A it A{, etc. 
are definite functions of m of the order m? at least. Their values are 

1 91 

A=-^-{4H+ (0}, ^ = -^-^p] + 3(i)}. 

^, i^', (7,, D' t , F'i, G{, Hi are got by merely interchanging [*] and (i) in the 
expressions for the corresponding undashed letters. Hence 

A _ 5m 2 2i 2 — 10a — 6 — 4 (t + 6) m — 27m 2 
i_ 16a 2 * 2(i 2 — 1) — 4m-f m 2 

^, _ _ 5m 2 14a 2 — 22i + 6 — 4 (7t — 6) ffl + 27m 2 

* 16**" 2(i 2 — 1) — im + m 2 

B = — — 6i 8 — 14a — 2 — 4 (3i + 2) m — 9m a 

* 16i 2 ' 2 (i 2 — 1) — Am + to 2 

fi'- _ 3ffl 8 10a 2 — 18t + 2 — 4 (5t — 2) + 9to 2 

* 16i 2- 2(t» — 1) — 4m + wi 2 
etc. 

Also we have 

Zl= i~Hr i(A + (^+-i-{Wi + (^M] + ••• •• 

The coefficients of £* in m 3 , s s , etc., are now to be obtained. From the forms 
of u, s we have, if the coefficient of £* in tt p s 4 be denoted by (u v s q ) t as before, 

(m p s ? ) < = 22 . . . . a/i^a; • . . . (p factors) X a„a t a w . . . . (q factors) , 

for all integral values, positive and negative, of j, k, J . . . . s, t, w . . . . con- 
sistent with the condition 

(j> + j + k + I ....) — (q + s + t + w ....) = »• 
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From this equation we see that 

(V?ffl) t = («*»«)_«• 

The coefficient of i? is more easily obtained by taking the coefficient of £ A , where 

h = i— p + q, 

and therefore 

(j + k + I +....)_ ( a + « + w + ....) = h. 

III. 

Tfte Coefficients « x an^ a_j. 

We have now obtained all the expressions necessary for the determination 
of the coefficients a t . But before developing them, some remarks must be made 
on the equations for o x and a_ x . 

It will be noticed that every term in equation (5) except the principal one 
is divided by the expression 2 (i 3 — 1) — Am + m 2 . When i = ± 1 , this reduces 
to — 4m + ot 2 , and the order of a x and a_ x is thus lowered by one power of m . 
The consequence of this is a great increase in the difficulty of obtaining their 
expressions. We have to carry them one order higher to get the same degree of 
approximation, and when we obtain them by the method of successive approxi- 
mation, each process, instead of carrying our expressions two orders higher, takes 
them only one ; the number of such processes is consequently doubled^ and the 
numerical multipliers of the various powers of m are very much increased in 
complexity. 

A second disadvantage arises from another cause. If we form the expan- 
sions of a 1 and a_ x in ascending powers of m from equation (5), we find large 
and continually increasing multipliers of the various powers of m which tend to 
lessen the value of the successive approximations when expressed numerically, 
to such an extent that the series, instead of proceeding in a progression whose 
ratio is roughly m or 1/12, proceeds in a progression with a ratio of about 5m 
or 2/5. As this ratio seems somewhat regular between the successive powers, 
I have used the following method to discover and counteract the effect of the 
slow convergence. In the following, squares and higher powers of l/a! are 
neglected. 
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Going back to the equations (4) obtained directly from the equations of 
motion, put **= 1 in them, and write down the coefficients of a x , a_ 1 , neglecting 
terms of the seventh and higher orders. These principal terms then become 

a^ [3 + 6m 4- \ m 2 4- (? + 10m + ~ m 2 ) -|- + \ m 2 . ^] 

+ a Q a_ x [l 4- 2m 4- 9m 2 4- (1 — 2m 4- 9m 8 ) ~] , 

a oai I— 3 — 2m — (5 + 2m) -?L + 3m 2 ^=*1 

L a a J 

+ a a_i ["— 1 — 2m 4- 3m 2 + (1 — 2m + 3m s ) ^=il . 

Now a 2 /a and a_ 2 /«o involve m only, and their expressions have been found 
by Mr. Hill. He gives 

^2 3 „ 1 „ 7 . 

— = _„,■+_„• + _.„»« + 

«- 2 19 5 , 43 . , 

Substituting these values and expanding, we have for the principal terms, 
«o«i [ 3 + 6m + -jg- m 2 4- ] 4~ «o«-i [ 1 + 2m + -— m 2 -f ], 

a Q a x [— 3 — 2m — — m 2 + ] 4- « «-i [— 1 — 2m + -^- m 3 4- 1 . 

Denote these expressions by a^P + a ( /i_ 1 P l and a^Q -+- aociLjQ' respectively. 

When we wish to separate out % and a_ x by the ordinary process of solu- 
tion, i. e. multiplying the equations by Q', P' and subtracting, also by Q, P and 
subtracting ; the principal terms become 

a <h(Pg-P'Q), 
a a^(QP-Q'P). 

The resulting expressions for a x and a_ x thus found would therefore be both 
divided by the factor PQ 1 — Q'P. Expanding it in powers of m, we find these 

expressions will be multiplied by 

— [ 1 — 4m m i — .... I . - — , 

\ 8 /4m 

or b J 1 f ■> , , ■ 165 2 , \ 
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We thus have an explanation of the slow convergence of the series, and by 
retaining this divisor, we shall be able to obtain series which will converge much 
more quickly. By proceeding in the following way, the values of a x and a_ x 
can be formed to the seventh order in one approximation ; and it will also be 
shown that to this order, our values will be sufficient to obtain the coefficient of 
this part of the Parallactic Inequality from the expansions, with an error less 
than two hundredths of a second of arc. 
In equations (4) put i = 1 , we have 

2/¥V-i [f +i + 1 + 2m (2/+ 1) + -§- m*] 

9 
H-S^m* [a j _ 1 a_ } + a_ ) _ l a j _^ — — r^L^ 

g 
Sja/f/-! [— 2/ — 1 — 2m] + 2, -g- m 2 [a j _ 1 a_ j —a„j_ 1 a j _ 2 ~] = — wiWj. 

Multiply the second equation by 2m , add to the first and write out the result- 
ing equation and the second equation, for errors of the sixth order at most in m, 

^['+4 + t( 7+ ^) + ^(i rf + 6 *" , ) 

+ °W ( 7 + "x)] = ~ m * (Zl + 2mMl)l (6) 

a^ [3 + 2m + -^- (5 + 2m) — 3m 8 ^ + 3m 2 ^1 

+ «o«3 [^ (7 + 2m) + ^ (9 + 2m)] 

+a a_ 1 [l + 2m — 3m 2 + ^ (— 1 + 2m + 3m 2 )] 

+a a_ 3 [3m 8 - 3m 2 -J- + (- 3 + 2m) «=1 + (- 5 + 2m) ^=- 4 ] = rfl, 

The equations (6) will be used for a x and a_j instead of those obtained from 
the general form (5). 
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IV. 

The Determination of the Parts of a { depending on the First Poicer of the Ratio 

of the Mean Distances. 

The coefficients a 4 will now be obtained to the order »i 5 .a /a'. To obtain 
them to this order we have 

(« 3 )i = 3a 2 a_ 2 , (u 2 s\ = a% + 2o (af + aL 2 + a z a_ 2 ) , 

(s 3 )! = 3ata_ i + 3a!_ 2 cr , («« 2 )i = «o (2a_ 2 + a,) , 

(u 3 ) 3 = (s 3 )_ 3 = a 3 ,, (« 2 «) 3 = (ws 3 )_ 3 = og (2a 2 + a_ 2 ) , 
(w 3 ) 5 =(s 3 )_ 5 = 3aga 2 , 

the rest of the coefficients of £ in it 3 , s 3 , u~s, s"u being zero to the order taken. 
Whence we have 

L\ = ~^r [-y («o«-s+«oa-4+«-2) + "y («o+2a 2 +2ai 2 +2a 2 a_ 2 +2a a.. 2 +cr a 2 )] , 
-M"i= -^r [-g- (a «-2— «o«-4— «iz) + y («o+2al+2a:L 2 +2a 2 a_ 2 — 2a a_ 2 — a,^)] 

for substitution in the equations (6). Also the right-hand side of equation (5) 
becomes when 

i= 3, -£- [A s a% + B 3 a (2a 2 + a_ 2 )], 









3 , -4- [ALA + BL 3 a (2a 2 + «_,)] , 



a, 



5 , —j- A$ . 3a a 2 , 



a 



t = — 5, ~7-^ 5 .3a a 2 , 



a 



whence, writing out equations (5) for a s , a_ 3 , a 5 , a_ 5 with a sufficient number of 
terms to obtain these quantities as far as the order m 5 . a /a', we have 



«o«3 = [3, 1] a!a_ 2 +[3, 2] a_ t a 2 + [3] 2a 1 a + — r [^s«o+ 5 3«o(2a 2 +a- 2 )] 

+ [3,4] a^ + [3, — 1] a_ja_ t + [3] 2a 2 a_ 1( 
a «-3 = [— 3, — 2] aja_ 2 + [—3, — 1] a_ 1 a i + (— 3) 2a 1 a 

+ -^rlA'_ 3 at + BL 3 a (2a, + a_ 2 )-] 

+ [— 3, 1] a^ + [— 3, — 4] a^a^i + ( — 3) 2a 2 a_ 1) J 



(7) 
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a a 5 = [5, l]a 1 a_ i + [5, 4] a i a_ 1 + [5, 2] a % a_ 3 

+ [5,3] a<p,^ 4- [5](2a 3 a + 2a a « 1 ) + -£- J. 5 .3a a 3 , 
a a_ 5 = [— 5 , — 4] a 1 a_ i + [— 5 , — l] a k a_ x 4- [— 5 ,. — 3] a % a_ z 

ft 

+ [—5, — 2]a 3 a_ 2 +(— 5)(2a 3 a +2a iJ a 1 ) + —7- ^1 5 . 3a a 2 . 



K8) 



The usual method is to solve the equations (6) for a first approximation of a lt a_ x 
by neglecting d s , a_ 3 , obtain a s , a_ 3 from (7) with these values of a lt a_ x and pro- 
ceed then for a second approximation of a lt a_ x . I prefer to do this in one pro- 
cess as follows. 

Since a 2 , a_ 3 , a 4 , a_ 4 are known expansions of m in series, we can to the 
order given, write the equations for a 3 , a_ 3 in the form 



a s = aa! 4- /?a_i + y -A- a , | 



an 



> (9) 



a_ 3 = a'rtx 4- /3'a_! + / — 7 



a' 



«o. 



where a, @ , y, a', @', y' are known functions of m. Substituting these values of 
a 3i «-3 in equations (6) for a b a_ x , we obtain linear equations of the form 



a x .x 4-a_ x .^ =^.-^r«o. 
aj.x' +a_ 1 .^' =(i'.-4- a , 

(J/ 



(10) 



where «, a, ^, x', Jl/, ^' are known functions of m. This method of procedure 
has the advantage of giving us the whole value of the functions of m used when 
dealing numerically with the equations; while in finding their expansions alge- 
braically, we lose much less than by the usual method. Also, if we wish to pro- 
ceed to a higher approximation, by using the equations always in the forms (9) 
and (10), we shall merely get terms in (9) independent of a x , a_i, and conse- 
quently in (10) terms added only to the right-hand sides of the equations. Then 
the coefficients of a lt a_ x in (10) will remain unchanged, and the small changes in 
(i, (i' can be easily obtained. This will not be found necessary here, as the 
results obtained up to the order indicated above are sufficiently accurate. 



y(ii) 
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Expressed in the form (9), equations (7) become 

,/51 ,3371 13451 \ ,/15 ,167 607 A 

, ,/ 5 , 45 .' 645 1 781 A a o 
+ »»*(-gr + -gT™ + gTsF*" 312 s m 7V "<" 

„/ 25 , 103 , 2659 A ,/ 3 19 215 A 

a_ 3 = m 8 (^— ^r + ^-^ m + gg-^io wi 8 J a x — W^-gf + "gT m + 3^ TO j «-: 

/55 175 229 , 8455 A a 

-™^-gr +- 2 T- OT - 3 T2 riOT 312 s m j¥" a « 

and therefore from equation (6) we can obtain 
/ 29 7 3 1163 4 209893 \ 

^V +16" + T O 10 m ¥^W ) 

/ 61 a , 13 , 68563 , 338477 A 

+ a ->( 1 +U m +-3~ m ~^2 i5m 3^2^*"; 

a / 3 , 3 , 12795 . 96643 . 752323 A 
= — -£- • «o ^-g- m + -4- m $tr- n»* j^ - m O 1 ^ m J ' 

/o ■ „ . 15 , , 23 , , 20645 . , 362467 A 
a x {% + 2m + -jg- «n* + -y m 3 + -g-^io ™ 4 + ^ gl> m 5 J 

/ 29 , 17 , 55379 , 891109 .\ 
+ a_ 1 (^l + 2m — ^g-m 3 — -gj-Tn 3 — -g^-gjo m 4 3372^ m J 

a„ / 3 „ 3165 . 21333 . 1300401 „\ 



(llo) 



Subtract the first of these equations from the second, the resulting equation is 
then divisible by m; also, multiply the second equation by 3 and subtract from 
the first. We have then two equations from which a lt a_ x must be found. 
Performing these processes, we get 



3 
a. 



! a 1 r 15 , 15 , . 123 . , 13975 . , 1709047 el *) 
- = ^-^FL-32 m ~16 m 3 + lu m +-W- m + ~0^ m J'l 

a_i « 1 T , 45 a , 21 3 2763 4 13449 5 2979511 -, [ ^ 
^ = ^-^L+32™ +16™ W m 2T m 0^ m J'J 



where mt, the common divisor, is given by 



37 3 17 89963 

mt—m — Am? — m 3 ■ g- m 4 — 38 2l0 w 5 . 



16 
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We can of course divide out by m. It has been kept in here to show in what 
manner a transformation of the form 



m = 



1 + am' ' 

if indicated by theory, must be made.* We should get different expressions for 
the coefficients a x and a_ x , kept still in the form (12) according as the factor m 
had been divided out or not. 

Using the values (12), after dividing out by the factor m, since it can now 
be any time replaced if necessary, equations (11) give 

<H l a T 5 . 45 . 545 2 781 3 

vn /45 , "8165 66251 .\-i 

- t \ ~r + o" u m ~ -$rw m ) J ' 

a_ 3 . a r 55 175 ,229 „ , 8455 . 

, m /15 . 551 7459 ,\n 

+ -TVW + s^ m — &&**)}' 

and thence from equations (8), 

«s «o ^rlOS , 797 2655 m -\ 

This method has the advantage also of giving numerical results along with the 
algebraical expansions, and thus we can obtain some idea of the errors which are 
produced by the portions of the series neglected when we expand the various 
functions used, in powers of m. In working out the numerical results by this 
method, it appeared that the principal parts of the errors in a x , a_ x were due to 
the neglect of the higher powers of m in solving the linear equations for a x and 
a_u and that even this part was not very great. For example, the value of t 

found was 

37 , 17 , 89963 , 
1 — 4m g- m a — m 3 — -grgnr »**. 

giving T=. 6443757. 

The value found by the numerical process to the same order is 

* = . 6444540, 
a difference of .0000783. 

*See Monthly Not. E. A. S., Vol. Ill, No. 2 . 
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The term 89963 



3 2 . 2 10 



m 4 = —.0004171, 



so that the neglected portion is less than one-fifth of the last term in t calcu- 
lated. 

The expressions given above, when transformed, agree with those given by 
Delaunay as far as the order rrt.a^d. The multiplier of in 5 .a /a' differs in all 
these coefficients by a small amount, causing, however, only a small difference 
in the numerical values of the coefficients of a few ten-thousandths of a second. 

V. 

The Determination of the Portions of the Coefficients depending on Powers of 

1/a' higher than the first. 

If now we wish to obtain the portions of the coefficients depending on 
(a /a'f and (a fa'f, very little extra labor is necessary. In the coefficients with 
even suffixes only even powers, and in those with odd suffixes, only odd powers 
of this ratio occur. Let $a r denote the portion to be added to a r depending on 
the square or cube of this ratio according as r is even or odd. The even coeffi- 
cients will be obtained below as far as the order m? (a /a'Y, and the odd ones as 
far as the order rn % (a /a') 3 . As Ba is of the order rn % (a /a'y, we have to the 
required degree of accuracy from equation (5), after obtaining the necessary 
terms in L, M, 

a fe 2 = [2, 1] a x a_ x + [2, 3] a 3 a x + [2, — 1] a_ x a_ 3 

+ ■%■ [3^a a_ 1 +5 2 a (2a 1 + a_ 1 )] + (jf) Vo, 
a<$a_ % = [— 2, — 1] a 1 a_ 1 + [— 2, l] a 3 a, + [— 2, — 3] a_ x a_ 3 



a 1 

2 



a. 



aMi = [ 4 > 3] asOL.! + [4, l]a 1 a_ 3 + -^-. ZA^a, + (-^hj #4«o> 
,3a_ 4 =[— 4, — ljaga^+f— 4, — 3] a&_ 3 +^ .ZAL&<fh+(^f) C-A, 
a 8a 3 = (-^r) # 8 «o, a 8a_ 3 = (jjfj G'_ 3 a%, 



(13) 
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For a 1 and a_ : we must use the equations (6). Differentiating these and 
remembering the orders to which our quantities are carried, we get 

a^a-^P + a 8a_ 1 P'+ 1a x ha % + a_^a_ z = — m*(8In -f 2tn8M 1 ), 

a$aiQ + a §a_iQ/-\- ba^ai — a_ x ha^ % = rrfhMx, 
where P, Q, P', QJ are the same functions of m which we had previously as 
multipliers of a t and a_j in their respective equations. Treat these in the same 
manner as we did the corresponding ones for a i} a_ 1? i. e. subtract the first equa- 
tion from the second and also multiply the second equation by 3 and add to the 
first, we get 

«o&*i(6 — P)+aM-i(Qf — P')=ni s l$L 1 +&Mi(l + 2m)]+ 2a 1 5a 2 +2a_ 1 5a_ 2 , 
a Sa 1 (SQ+ P)+a a 8a-i (3 Q+P) = m z [— *A+ &«i (3— 2rn)']—22a 1 $a i + 2a_ x ha~% , 
corresponding to the equations (11a), the coefficients on the left-hand side being 
the same. Also 

25 ,_ . s 45 _, , ,45 



dij = a (-^rj [-jg- (3a_! + a x ) + -^ 2 (a_ x + Oj) + -j^ «oj 
&Mi = a (~?-) [— (3a_! + oO + -gj a ] • 



Working out these and equations (13), and retaining the factor — wherever it 
occurs, we obtain the following series of values : 

fa, 3 /a \ 2 r 5 5 45.67 m 45 1 -| 

5a_ 2 . / a Vr 45 5 15.1253 m 45 .„ , „ . 1 1 

_ = m ^_ r j|__-^_ 0§m ^.__ + _ ir(5 _ 1 _ 3m) __J ) 

Sa 4 , / ff \ 2 r 7 , 7 15.25 m I 

£a_ 4 «/«oVr 7.17 21 , 15.45 ml 

&*i /«o\ 3 r 105 210 „ 1125 m 2 -\ 1 

5a_ x / a \3r315 ,270 , , 3375 m % -\ 1 

<Hi o / an \ 8 35 



m' 



■ffl- 



a -'" -VaV-3.2 11 ' 
£a_ 3 g / a \ V 35.43\ 

£ = -■(*)•!£.■ 
% , =--(^)'(-J^)- 



Motion which is a Function of the Mean Motions of the Sun and Moon. 157 



When these are compared with Delaunay's expressions they all agree as far as 
the order to which he has carried them. 

The results obtained by using numerical values from the outset are : 



= — .0641700 -^.-.086 (-%)' 
a' \ a' J 

^=+.1789909^.240 (£)' 

-2*- = _ .0000589 \ + .0001 (■%) 
a a' ' \ a / 

^=i = - .0029370 -%• - .0016 (^-Y 
a a' \ a ' / 

— = + .00000047 ■% + .0000 f -^-V 

^ = - .00000184 -^ - .0010 C-^-V 
^= + .0072*4^)' 
£ = +. 000000 (^)', 



a 



= -. 007482 (-£-) 
= -. 001428 (£.)' 



tion 



The quantity fe must now be obtained. For this purpose I use the equa- 

(? + 2mD - (=?)«=- 4 ■><•+•)-*. ^ • f 1 • 

Substituting in this 

and taking out the coefficient of £ we get with the former notation, 

* (&p)r ( 1+2wi + 4 ™ 2 ) «•+ 4- ™ 2 «-* + ^r [it <n + 4 ( %2 + 2MS )J 

+ "S 5 " ["16 (s3)l + "16" ^ + ~W (Ms2)l + T6~ ( " 2s)l ] + ' ' ' ' 
When the parallactic terms are neglected Hill finds that 



where 



r«7Yl + «i) a -i* 
a = a ]_ g ] ' 
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and a ff is the value of the right-hand side of the equation. Let Sa , SJ, S3 
be the corresponding quantities due to the parallactic terms. Then since 

J+U, H+m 
are defined in the same way as J, H, we have 

a + Sa -a [ E+ ^ H J • 

Neglecting powers of —f- above the second, we get 

a — 3 ' J 3 ' H ' 

We shall obtain Sa to the order m 3 . (a /a'f. For this 

a' L 8 a ' 4 a J a' 16 

r 9 225 m'-]/a \* 

= l-U m+ ^'-rXaT) • 

a % u (Xa P) a 2 

J+ 5/ is the coefficient of £ in 7—1 or of £° in 7^ — fS lu ; whence we get 



3 af-f aLi _9_ (q 1 + a_ 1 ) 3 

" 2~ ' <? ~ ~ " ' ~ ~^~~" 

m 2 675 



dJ - ~ 2 * a 2 + 4 * a 2 



T* 



(i + *-)(t-)- 



Also H=l + 2m+±-<* + ±nf^', 

21 
/= 1 + -gr »J 4 + 

/ rn % \ 

a = a{l g J , 

we therefore have 

5a / a V '(225 1 + 4m , 75 m . 3 .1 

and numerically, 

^=-^Y.00965. 
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VI. 

Transformation to Polar Coordinates. 

The inequalities will now be expressed in polar coordinates. For this pur- 
pose we have, if V be the true longitude of the moon, 

rcos(V-nt) = —Xa,^ + ?-«) 



= a (l+-i-2A^), 



rsin(7- nt) = ^—Xa,^-^*) 

1 



= ««^^T 2 ^' 



where -4 i .a = a 4 + a_ t and 5,-.a =6i — &_«, and consequently .i^J.^, 
JSi = — B_ i . And since 

6 s 
tan = - — \- . . . . , 

F -^= 2 v^iPAr][i-4- s ^+x^^ 2 -'--'] 

+ 247^1 [^T[l---..] 3 +.-.- 



The coefficient of sin D or (£* — ^ _1 )/2-</— 1 is therefore 

JBi + 4" 0M.-A*.) + 4" 0M.-4A) + 4"^(^- B ® + ■ • • • > 

or a x — a_ x a^a^^ — a^a_i + a i a_ 3 — a 3 a_ a , 2 (ax — «_ i) a^X- % 

a a uo 

and similarly for the other coefficients. Performing these operations we get the 

following expressions, omitting those terms dependent on rn only. 

In longitude, 

15 , 9 , 1951 , 41585 , 2096751 
■ xW + _^__ w t_.-_ OT 4_- i§nr , w , 35 4 _2345_ 

~ 37 I 17 I 89963 ~ ! 1024 W 12288 

1 — 4m — — r— mr ^— m 6 — nc>1/> nr ( = <r) -o".ooo6 -cr.0002 

8 6 921b v ' 



- 128".0545 



a' 2 1 64 ' T + 8 t + 128 ' T 3 j 

— C.0006 — 0".0001 — O".00O4 



a . n 
— rsin D 
a' 



160 Brown : On the Part of the Parallactic Inequalities in the Moon's, etc. 

, l~5»i 2 /c , „ n . 225 o(5 + li™) , 15.747 m s -1 a 2 . on 
+ L^(5+8 W )+ I ^^L_^ + __._ r J_ rSm2Z > 

C.0078 +0".0276 +0"0060 

, r 15 2 , 55 3 41 , 2309 . 1 /255 , , 7543 , 134525 _\ 
+ L-32 m + 128 m - 1536 m ~ 9216 m ~ T (l28 m + 1536 W 8192" m ) 

1".6172 +0".1199 -0".0006 -0".0005 — 0".9872 

. a 2 385 .-1 a . 
a' 1536 J a' 

+ 0".0000 

, r63 „ . 77 „ 1725 m s -t a 2 . „ 
+ 1.256 m + 480 m ~ 2048 * TJ ^ sm4Z) 

0".0021 +O".0O01 — 0".0008 

, r75 4 2797 5 18415 m* a 3 777 ,-1 a . B _ 
+ Ll28 m *+ 2048 m - -8191- ' — + ^ ' 5120 W J ^ Sm 6i? ' 

0".0132 + 0".0025 — 0".0064 +0".0000 

In parallax, 

1 r / 15 3 . 1869\ 1 -1 a _ 

^L-(,T6- m + -8" m -^56-;vJw- cosZ? 



, 1 r 25 65 , 459 m 3 -i 1 n _ 
+ -^L-64 m + 256 «'- 256 ' T"J ^ C ° S ^ 

0".0223 +C.0O12 — 0".0139 

None of the other parts produce in parallax, coefficients so great as one thou- 
sandth of a second. 

The numerical values of these coefficients are obtained by using 

m = .0808489338, 
a 



, —.00255878, 
a' 



If we use the numerical value of m from the outset, we obtain the more 
accurate series of coefficients : 
In longitude, 

— 128".070 sinD+0".039 sin 2D+0".750 sin 3D+ 0".001 sin4Z>+0".008 sin 52), 

and in parallax 

— 1".001 cos D + 0".008 cos 32). 

For the discussion and comparison of these results with those obtained by 
Delaunay, see Monthly Notices of Royal Astron. Soc. f Vol. LII, No. 2. 
November, 1891. 



